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Abstract Recently, using Greenwich and Solar Optical Observing Network 
sunspot group data during the period 1874-2006, Javaraiah (Mon. Not. Roy. 
Astron. Soc. 377, L34, 2007: Paper I), has found that: (1) the sum of the areas 
of the sunspot groups in 0° - 10° latitude interval of the Sun's northern hemi- 
sphere and in the time-interval of —1.35 year to +2.15 year from the time of the 
preceding minimum of a solar cycle n correlates well (corr. coeff. r = 0.947) with 
the amplitude (maximum of the smoothed monthly sunspot number) of the next 
cycle n+1. (2) The sum of the areas of the spot groups in 0° - 10° latitude interval 
of the southern hemisphere and in the time-interval of 1.0 year to 1.75 year just 
after the time of the maximum of the cycle n correlates very well (r = 0.966) 
with the amplitude of cycle n + 1. Using these relations, (1) and (2), the values 
112 ± 13 and 74± 10, respectively, were predicted in Paper I for the amplitude of 
the upcoming cycle 24. Here we found that the north-south asymmetries in the 
aforementioned area sums have a strong w 44-year periodicity and from this we 
can infer that the upcoming cycle 24 will be weaker than cycle 23. In case of (1), 
the north-south asymmetry in the area sum of a cycle n also has a relationship, 
say (3), with the amplitude of cycle n+1, which is similar to (1) but more 
statistically significant (r = 0.968) like (2). By using (3) it is possible to predict 
the amplitude of a cycle with a better accuracy by about 13 years in advance, and 
we get 103 ± 10 for the amplitude of the upcoming cycle 24. However, we found 
a similar but a more statistically significant (r = 0.983) relationship, say (4), by 
using the sum of the area sum used in (2) and the north-south difference used 
in (3). By using (4) it is possible to predict the amplitude of a cycle by about 
9 years in advance with a high accuracy and we get 87 + 7 for the amplitude of 
cycle 24, which is about 28% less than the amplitude of cycle 23. Our results also 
indicate that cycle 25 will be stronger than cycle 24. The variations in the mean 
meridional motions of the spot groups during odd and even numbered cycles 
suggest that the solar meridional flows may transport magnetic flux across the 
solar equator and potentially responsible for all the above relationships. 
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1. Introduction 

Prediction of the strength of a solar cycle well in advance is important for 
predicting the space weather because solar activity affects space weather in 
several ways (Hathaway and Wilson, 2004; Hathaway and Wilson, 2006; Kane, 
2007a). In addition, it may help for understanding the basic physical processes 
responsible for solar cycle (Dikpati, de Toma, and Gilman, 2006; Choudhuri, 
Chatterjee, and Jiang, 2007; Cameron and Schussler, 2007; Jiang, Chatterjee, 
and Choudhuri, 2007). Many characteristics of a solar cycle and of its previous 
cycles have been used to predict the amplitudes of the cycle (Hathaway, Wilson, 
and Reichmann, 1999; Li, Yun, and Gu, 2001; Kane, 2007a; Obridko and Shelt- 
ing, 2008). The existence of the north-south asymmetry in solar activity is well 
known since several decades and has been extensively studied using the data 
on almost all the solar activity phenomena including sunspot number, sunspot 
groups, solar flares, prominences/filaments, photosphcric magnetic fields, so- 
lar rotation and differential rotation, meridional flow, etc., and revealed many 
characteristics of it (Roy, 1977; Swinson, Koyama, and Saito, 1986; Antonucci, 
Hocksema, and Scherrer, 1990; Garcia, 1990; Hathaway and Wilson, 1990; Yi, 
1992; Carbonell, Oliver, and Ballester, 1993; Varma, 1993; Sokoloff and Nesme- 
Ribes, 1994; Atac, and Ozgiic, 1996; Javaraiah and Gokhale, 1997; Duchlev 
and Dermendjiev, 1996; Li et ai, 2002; Javaraiah, 2003; Georgieva and Kirov, 
2003; Georgieva et ai, 2005; Ballester, Oliver, and Carbonell, 2005; Gigolashvili, 
Japaridzc, and Kukhianidzc, 2005; Knaack, Stenflo, and Berdyugina, 2005; Joshi 
and Pant, 2005; Temmer et ai, 2006; Zaatri et ai, 2006; Javaraiah and Ulrich, 
2006; Chang, 2008). However, so far no physical or even a significant statistical 
relationship between the solar activity and its north-south asymmetry is known. 
Recently, Javaraiah (2007, hereafter Paper I) has found the following statisti- 
cally high significant relationships between the sums of the areas of the sunspot 
groups, An and As (normalized by 1000 ), in 0°-10° latitude interval of the 
Sun's northern hemisphere and in the time interval of : T m +(— 1.35 to +2.15), 
i.e. —1.35 year to +2.15 year from the time of the preceding minimum (T m ) - 
and in the same latitude interval of the southern hemisphere but in the time 
interval of T M : Tm + (1.0 to 1.75), i.e. 1.0 year to 1.75 year just after the time of 
the maximum (Tm)- during a solar cycle n and the amplitude (Rm, maximum 
of the smoothed monthly sunspot number) of the next cycle n + 1 : 

Rm.u+i = (1-72 ± 0.19)A N , n (T* ) + (74.0 ± 7.0), (1) 

Rm,u+i = (1-55 ± 0.14)A S ,„(T M ) + (21.8 ± 9.6), (2) 

where n = 12, 13, ,23 represents the Waldmeier solar cycle number. The 

corresponding correlation coefficients of these relations are r = 0.947 and r = 
0.966, respectively. By using Equations (1) and (2) in Paper I we have found 
112 + 13 and 74 ± 10, respectively, for Rm of the upcoming solar cycle 24. Since 
Equation (2) is statistically more significant than Equation (1), hence, in Paper I 
we have predicted the value 74 ± 10 for Rm of the upcoming cycle 24. In the 
present paper we have made a detailed study on the north-south asymmetries in 
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the sums of the areas of the spot groups, i.e., differences between the An,™ and 
As,n, during both the time intervals and T M of the solar cycles n = 12 to 
23. This enabled us to find two more new relationships, (3) and (4), and using 
these to improve the prediction in Paper I for the amplitude of solar cycle 24. 
The physical significance of all the aforementioned relationships is discussed. 

It should be noted here that in most of the studies of the north-south asymme- 
try of a solar activity phenomenon, the north-south asymmetry is determined 
using the conventional formula, , where N and S are the corresponding 
quantities of the activity phenomenon in the northern and the southern hemi- 
spheres, respectively. This fraction has a strong 11-12 year periodicity (Car- 
bonell, Oliver, and Ballester, 1993). Yi (1992) pointed out that in the north-south 
asymmetry, derived by using this formula, the 11-12 year periodicity has no 
statistical significance because mostly it is an artifact of the 11-year cycle of 
N + S. In order to verify this Javaraiah and Gokhale (1997) determined the 
power spectra of both the N — S and the of the sunspot number data 

and found that a peak at 11-12 year in the power spectrum of N — S is 
statistically insignificant, whereas it is statistically very significant in the spec- 
trum of jlfT§ i confirming the doubt of Yi (1992). Recently, Ballester, Oliver, 
and Carbonell (2005) also confirmed the same. Therefore, the difference N — S 
seems to represent the north-south asymmetry of a solar activity phenomena 
more appropriately than the fraction j^lf ■ In addition, in the present analysis 
we found that the correlation (r = 0.90) between of a cycle n and the 

Rm of cycle n + 1 is much weaker than that (r = 0.968) between the difference 
An — As during T£ of a cycle n and the Rm of the cycle n + 1. Therefore, here 
we have used the difference An — • 

Since a large number of abbreviations are used here, hence for the sake of 
the readers convenience we listed below the meanings of all the abbreviations 
(thanks are due to the anonymous referee's suggestion): 

• n - the Waldmeicr solar cycle number, 

• T m - the preceding minimum epoch of a solar cycle, 

• Tm - the maximum epoch of a solar cycle, 

• R m - the value of smoothed monthly sunspot number in T m , 

• Rm - the value of smoothed monthly number in Tm, 

• Tj* - the time-interval of —1.35 year to +2.15 year from T m , 

• T M - the time-interval of 1.0 year to 1.75 year just after Tm, 

• AN,n (T^) - the sum of the areas of spot groups in 0° — 10° latitude interval 
of the northern hemisphere during Tjjj of a cycle n, 

• As tTl (T^) - the sum of the areas of spot groups in 0° — 10° latitude interval 
of the southern hemisphere during Tj£ of a cycle n, 

• ^N,n(T M ) - the sum of the areas of spot groups in 0° — 10° latitude interval 
of the northern hemisphere during T M of a cycle n, 

• As, n (T M ) - the sum of the areas of spot groups in 0° — 10° latitude interval 
of the southern hemisphere during T M of a cycle n, 

• SA n - the difference A N ,„(T*) - As,JT£), 

• ^NS,n - the sum A Nt „(T£J + A Sj „(T M ), and 

• A SA s,n - the sum SA n (T*J + A s , n (T M ). 



javaraiah_msrev.tex; 8/03/2009; 7:44; p. 3 



J. JAVARAIAH 



In the next section we describe the data analysis and the results. In Sec- 
tion 3 we discuss about the implications of the relationships among A-^ t „ (T,* ), 
As, n (T M ), SA„(T^), A s asm, and Rm,u+i for understanding the underlying phys- 
ical process of the solar cycle. For this purpose we have determined the variations 
in the mean meridional motions of the spot groups in the northern and the 
southern hemispheres during the odd numbered cycles (11, 13, 15, 17, 19, 21, 
23) and the even numbered cycles (12, 14, 16, 18, 20, 22), using the superposed 
epoch analysis as described in Javaraiah and Ulrich (2006). In the same section 
we compare our prediction with other authors' predictions for the Rm of solar 
cycle 24, particularly with those predictions based on the flux-transport-dynamo 
models, and the spectral analyses and magnetic oscillations models. In Section 4, 
we have given the summary of the results and the conclusions. 

2. Data Analysis and Results 

2.1. Data 

Here we have used the values of and As, i.e., the sums of the areas of the 
sunspot groups in 0° — 10° latitude intervals of the northern and the south- 
ern hemispheres during the time intervals Tjjj and T M of cycles 12-23, which 
were determined in Paper I by using the Greenwich sunspot group data during 
the period 1874-1976, and the sunspot group data from the Solar Optical 
Observing Network (SOON) of the US Air Force/US National Oceanic and 
Atmospheric Administration during 1977 January 1-2005 September 30. The 
recently updated these data were taken from the NASA web-site of David Hath- 
away (http : / / solarscience .msf c .nasa.gov/greenwch. shtml). In Table 1 we 
have given the values of An and As during both the time intervals X 1 ,* and T M , 
and the maximum (Rm) and the minimum (i? m ) amplitudes (the largest and 
the smallest smoothed monthly mean sunspot numbers) and the corresponding 
epochs Tm and T m of the solar-cycles 12-23 (the values of T m , Tm, R m and Rm 
of cycles 12-23 are taken from the website, ftp://ftp.ngdc.noaa.gov/STP 
/ SOLARJDATA/ SUNSPOTJfUMBERS) . The details of the data reduction and the 
analysis were described in Paper I. The method of determination of An(T 1 *) 
and As{T M ), as described in Paper I, is as follows: first we determine the sums 
of the daily areas of the spot groups in the time-intervals whose lengths are 
chosen arbitrarily around the minimum and the maximum of a cycle n and Rm 
of the cycle n + 1 or n + 2, etc., and then we determine the exact time intervals 
by increasing or decreasing the arbitrary intervals with 0.05 year until we get 
the maximum correlation. We have chosen this method because the lengths of 
the solar cycles vary. However, it may be possible to obtain the similar results 
by binning the available data for the entire period (1974-2007) into successive 
intervals of suitable sizes and then picking the intervals close to the minimum 
and maximum epochs of the cycles. 

Note: In case of 0° — 10° latitude interval of the southern hemisphere, a 
reasonably good correlation is found between the sum of the areas of the spot 
groups during a long time-interval approximately one year after the maximum 



javaraiah_msrev.tex; 8/03/2009; 7:44; p. 4 



PREDICTION OF THE AMPLITUDE OF A SOLAR CYCLE 



Table 1. The maximum (-Rm) an d the minimum (R m ) amplitudes (the largest and the smallest smoothed monthlyl 
mean sunspot numbers) of the solar-cycles n = 12-23 and the sums of the areas of spot groups (normalized by 1000) inl 
0° — 10° latitude intervals of the northern hemisphere (An) and the southern hemisphere (As) during the time intervals 
= T m + (—1.35 to 2.15) and = Tm + (1.0 to 1.75), where Tm and T m represent the maximum and the preceding 
minimum epochs of the solar cycles, respectively. 



Cycle 


Minimum 


Maximum 




Around minimum 






Aftcr maximum 




n 


T m 




Tm 


Rm 


T* 

m 


A N 


A s 


T* 


A N 




12 


1878.9 


2.2 


1883.9 


74.6 


1877.55 


-1881.05 


9.47 


9.84 


1884.90- 


1885.65 


41.38 


42.11 


13 


1889.6 


5.0 


1894.1 


87.9 


1888.25 


-1891.75 


3.22 


30.83 


1895.10- 


1895.85 


18.67 


32.64 


14 


1901.7 


2.6 


1907.0 


64.2 


1900.35 


-1903.85 


12.98 


10.94 


1908.00- 


1908.75 


35.75 


54.64 


15 


1913.6 


1.5 


1917.6 


105.4 


1912.25 


-1915.75 


3.74 


7.93 


1918.60- 


1919.35 


84.76 


34.58 


16 


1923.6 


5.6 


1928.4 


78.1 


1922.25 


-1925.75 


33.96 


13.03 


1929.40- 


1930.15 


64.20 


75.96 


17 


1933.8 


3.1 


1937.4 


119.2 


1932.45 


-1935.95 


29.96 


7.43 


1938.40- 


1939.15 


49.12 


82.01 


18 


1944.2 


7.7 


1947.5 


151.8 


1942.85 


-1946.35 


69.35 


12.06 


1948.50- 


1949.25 


70.31 


119.65 


19 


1954.3 


3.4 


1957.9 


201.3 


1952.95 


-1956.45 


15.23 


13.18 


1958.90- 


1959.65 


116.12 


53.01 


20 


1964.9 


9.6 


1968.9 


110.6 


1963.55 


-1967.05 


50.31 


6.89 


1969.90- 


1970.65 


37.96 


78.28 


21 


1976.5 


12.2 


1979.9 


164.5 


1975.15 


-1978.65 


60.05 


25.06 


1980.90- 


1981.65 


58.42 


83.53 


22 


1986.8 


12.3 


1989.6 


158.5 


1985.45 


-1988.95 


29.85 


20.93 


1990.60- 


1991.35 


41.33 


67.48 


23 a 


1996.4 


8.0 


2000.3 


120.8 


1995.05 


-1998.55 


21.99 


20.17 


2001.30- 


2002.05 


77.27 


33.58 





a indicates the incompleteness of the current cycle 23. 



to approximately two year before the end of cycle n and Rm of cycle n + 1 . Since 
the correlation is maximum for T M , hence we have considered it. Although 
is much shorter than T,* , As(T M ) is considerably larger than A^(T^). Hence, in 
spite of T M is short the AsiT-^) time series is also well defined. 

2.2. Variations in the Sums of the Areas of the Spot Groups 

Figure 1 shows the variations in A N (T,*), A S (T* 1 ), A N (T M ) and A S (T^) during 
cycles 12 - 23. In this figure we have also shown the variation in Rm of cycles 12 - 
23. In Figure 1(a), it can be seen that there exist considerable differences between 
An(T^) and As(T^). The amplitude of the cycle-to-cycle variation of An(T^) is 
much larger than that of the corresponding variation of A$ (T^). There is a strong 
suggestion of the existence of a periodicity of about «44 years (double Hale 
cycle) in A^iT^). Tj[j of the largest amplitude cycle 19 is one of the minimum 
epochs of the «44 years cyclic variation. A^(T^) leads Rm by about 13 years 
[as already known from the relationship, Equation (1) above]. The cycle-to- 
cycle variation in As(T*J is weak [A$(T*A is almost constant from cycle 14 
to cycle 20], but there is an evidence for the existence of a weak 80-90 year 
cycle (Gleissberg cycle). In Figure 1(b) it can be seen that there is no much 
difference between the maximum amplitudes of the variations in AnC^m) ano - 
As(T M ). As we already know from Paper I the pattern of As(T^) is similar to 
that of A N (T^) (r = 0.94), and A S (T M ) leads R M by about 9-10 years. A S (T^) 
has a strong « 44- year periodicity. A-^^-^) varies approximately in phase with 
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Rm {r = 0.6) and has a somewhat different periodicity, 33-44 year (33-year 
periodicity associates with the high level of activity). The correlation between 
A^(T M ) and Rm indicates that there is an influence of the latter on the former. 
It should be noted here that Rm of a cycle n leads A-^(T M ) of the same cycle by 
1.0 year to 1.75 year. There is no significant correlation either between A^(T^J 
and A S (T^) or between A N (T^) and A S (T M ). 

2.3. Variations in the North-South Asymmetries of the Area Sums 

Figure 2(a) shows the variations in the north-south differences 6A(T*^) = A^(T^)— | 
A S {T M ) and 6A(T M ) = A N (T M ) - A S (T M ) during the cycles 12-23, and Fig- 
ure 2(b) shows the corresponding cross-correlation functions, CCF(Rm, SA), of 
the cross-correlations between SA and Rm (a positive value of lag indicates 
Rm leads 5 A). In Figure 2(a) we have also shown the variation in Rm during 
cycles 12-23. In this figure it can be seen that there exits an approximate 
anticorrclation (r = —0.57 ) between SA(T^) and SA(T M ). In the same figure it 
can also be seen that there is a strong suggestion of the existence of the «44- 
year periodicity in both SA{T^) and 5A(T M ). A 33-44 year periodicity is also 
found in the power spectrum analysis of the data on north-south asymmetries in 
sunspot activity (Javaraiah and Gokhale, 1997; Li et at, 2002; Ballestcr, Oliver, 
and Carbonell, 2005) and long-lived solar filaments (Duchlev and Dcrmendjiev, 
1996). This periodicity is dominant in the north-south asymmetries of the solar 
equatorial and the differential rotation rates determined from sunspot group 
data (Javaraiah and Gokhale, 1997), and it seems to be present prominently in 
the climatically related phenomena and in the Earth rotation rate, which may 
be related to the 44-year cycle (double Hale cycle) of solar magnetic activity 
(Fairbridge and Hillaire-Marcel, 1977; Georgieva, 2002). 

2.4. The new and the Improved Predictions 

In Figure 2(b) it can be seen that the cross-correlation function CCF(Rm,SA) 
of Rm and SA(T*A has a strong peak (0.9) at lag = —1, suggesting that 5A(T£) 
leads Rm by about 13 years, which is equal to the phase difference between 
An(T^) and Rm- Therefore, using 8A(T^) of cycle n we can predict Rm of cycle 
n+ 1. The CCF(i?M, SA) of Rm and SA(T M ) has a maximum negative value also 
at lag = —1, but its magnitude is inadequate for predicting Rm.u+i by using 
SA n (T M ). 

The correlation between 6A n (T^) and Rm.u+i is high (r = 0.968, correspond- 
ing confidence level is > 99.99) for cycles n = 12 to 23. The corresponding linear 
regression fit between 5A n (T^A and Rm.u+i is: 

#m,»+i = (1-65 ± 0.14)<L4„(T* ) + (99.8 ± 3.9), (3) 

where uncertainties in the coefficients are the formal la (standard deviation) 
errors from the linear least-square fits. 

Using Equation (3) the amplitudes of the upcoming sunspot cycles can be 
predicted by about 13 year advance. The result of the least-square fit is shown 
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16 18 20 

Solar cycle number 



24 



Figure 1. Plots of the sums of the areas of sunspot groups - Ajq (filled circle-solid curve) and 
As (open circle-dotted curve) in 0° — 10° latitude intervals of the northern and the southern 
hemispheres, respectively, during (a) and (b) - versus solar cycle number. In both (a) 
and (b) the cross-dashed curve represents the variation in Rm- 
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Figure 2. Plots of (a) the north-south differences (<5A), An — Ag for the values of An and 
As given in Table 1, versus cycle number, and (b) the corresponding CCF(i?M> SA) versus 
lag, during the solar cycles 12-13 (Note: a positive value of the lag indicates that Ryi leads 
the 5 A). The filled circle-solid curve and open circle-dotted curve represent the 8 A and the 
corresponding CCF(Rm , 5 A) during T t * and T^, respectively. In (a) the cross-dashed curve 
represents the variation in Rm- 
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in Figure 3. The correlation between simulated amplitude (Pm) and Rm, and 
also the level of the significance, are found to be the same as that between SA n 
and Rm.u+i- Using Equation (3) we obtained the value 103 ± 10 for Rm of the 
upcoming cycle 24. 

It should be noted that since Rm of cycle 23 is already known, hence it is 
included in the fittings of Equations (1) and (3), although it is not needed for 
predicting Rm of cycle 24. That is, using the values of An(T^) and SAtT^) 
of cycle 23 and the linear-relationships determined from the 9 pairs of data 
points correspond to the cycle pairs 12,13 to 21,22, we can get almost the same 
values for Rm of cycle 24 as found above. Since cycle 23 will end soon, using 
Equations (1) and (3) it will be possible to predict an approximate value for Rm 
of cycle 25 in about 3 years time. Rm of cycle 24 is not needed for this. 

Since variation in Ag(T,* ) is weak (approximately constant), hence the pattern 
of the variation of the corresponding 5A(T^) is almost the same as that of 
An(T,* ) (r = 0.95 between these quantities). Therefore, the value obtained using 
Equation (3) above for Rm of cycle 24 is approximately equal to the value 
obtained in Paper I by using Equation (1). However, the statistical significance 
of Equation (3) is better than that of Equation (1). 

All the values, 112 ± 13, 74 ± 10 and 103 ± 10, obtained for R M of cycle 24 
using Equation (1), (2) and (3), respectively, are considerably (7%, 39% and 
15%, respectively) less than the value of Rm of cycle 23 and consistent with the 
suggestion that the level of activity is now at the declining phase of the current 
Gleissberg cycle (Javaraiah, Bertello, and Ulrich, 2005b). 

The statistical significances of the value, 103 ± 10, predicted here by using 
Equation (3) and the value, 74 ± 10, predicted in Paper I by using Equation (2) 
are equal, but the difference between these values is considerably large (it is 
marginally less than 2a). Therefore, we are unable to predict that between these 
values which will become close to the real Rm of cycle 24, although we are 
confident about that the cycle 24 will be smaller than cycle 23. 

On the other hand, Equations (1), (2) and (3) may represent the different 
components of the total contributions of An.t^T^), SA^Tj^) and As :tl (T^) to 
Rm.u+i- Hence, it is necessary to find a relationship between the combinations 
of A N ,„(T n * ), 6A n (T£) and A s ^ n {T M ) and the R M ,n+i- We find that the Ans,n = 
^N.n^m) + As, n (T M ), i.e. sum of the area sums used in (1) and (2), of cycle 
n correlates very well (r = 0.973) with the i?M,n+i- Further, the correlation 
between A 5j4s ,n = SA n {T*A + A s , n (T M ) and R M , n +i is very high (r = 0.983). 
From this relationship we get the following relation: 

Ru,n+i = (0.83 ± 0.05)A MS ,« + (57 A ± 4.7), (4) 

whose statistical significance is very high, i.e., higher than those of all the rela- 
tions (l)-(3). Figure (4) shows a scatter plot of A^s,n and Rm,u+i- By using 
Equation (4) the prediction of the amplitude of a cycle is possible by about 9 
years in advance with a high accuracy. Using Equation (4) we get 87 ± 7 for Rm 
of the upcoming cycle 24, which is about 28% less than the Rm of cycle 23. 

Figures 5(a) and 5(b) show the plots of the simulated amplitudes (Pm), simu- 
lated by using Equations (1) - (4), and the corresponding differences, Rm — Pm, 
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Figure 3. Plot of the correlation between the north-south difference 8A(T^) = An — As 
during T£ of cycle n and _Rm of cycle n+1, where n = 12, 22 is the Waldmeier cycle number. 
Near each data point the corresponding value of n + 1 is shown. The solid line represents the 
corresponding linear relationship. The value of the correlation coefficient (r) is also shown. 
The values of An, Ag and Rm are the same which are given in Table 1. 



against solar cycle number 13-23. In Figure 5(a) we have also shown the varia- 
tion in Rm during these cycles, and all the predicted values for Rm of cycle 24. 
The uncertainties in the predicted values correspond to la values of Rm — Pm 
shown in Figure 5(b). In these figures it can be seen that the values of Pm and 
Rm closely agree each other, in the cases of Equations (2) - (3). This agreement 
is much higher in the case of Equation (4). 

2.5. Will Cycle 25 be Stronger Than Cycle 24 ? 

According to the well-known Gnevyshev and Ohl, or G-0, rule (Gnevyshev and 
Ohl, 1948) a preceding even numbered cycle is weaker than its following odd 
numbered cycle. However, occasionally the G-0 rule is violated. Figure 6 shows 
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the annual variation of sunspot number during 1610-2006. As can be seen in 
this figure, cycle 23 is an anomalous cycle in the sense that the cycle pair 22,23 
violated the G-0 rule. The cycle 5 is also an anomalous cycle (the cycle pair 
4,5 violated the G-0 rule) and it is followed by the week cycles 6 and 7. That 
is, a violation of the G-0 rule seems to be followed by a few small cycles (see 
also Javaraiah, 2005). Therefore, the violation of the G-0 rule by cycle pair 22,23 
also indicates that the next cycle 24 (probably even the cycle 25) will be a weak 
cycle. 

In principle the G-0 rule can be used for predicting Rm of an odd numbered 
cycle from that of its preceding even numbered cycle. But this is possible only 
when we know in advance that the even- and odd-numbered cycle pair will not 
violate the G-0 rule. For instance, in order to predict Rm of cycle 25 from that 
of cycle 24 (when it will be known) by applying the G-0 rule, it is necessary and 
essential to know in advance that cycle pair 24,25 will not violate the G-0 rule. 
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Figure 5. Plots of (a) the observed amplitudes (Pm) and the simulated amplitude (Pm) - 
simulated earlier using Equations (1), (2) and now using (3), (4) - and (b) the differences, 
Rm — flvii against cycle number. The cross represents the -Rm- The filled circle, open circle, 
square and diamond represent the Pm and the difference, Pm — ^m> corresponding to Equa- 
tions (1), 2, (3) and (4), respectively. In (a) the corresponding symbols at cycle 24 represent 
the values of Pm, simulated using Equations (1), (2), and (3), and represent the predicted 
amplitude of this upcoming cycle. 
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Recently, Javaraiah (2005) shown that it may be possible to forecast the epochs 
of the violations of the G-0 rule well in advance from the epochs of retrograde 
motion of the Sun about the solar system baryccnter. Except that so far no 
plausible method is available for predicting the violation of the G-0 rule. As can 
be seen in Figures 1(a) and 1(b) the values of A^(T^) and As(T M ) are higher in 
an even numbered cycle than in its following odd numbered cycle, say 'inverse 
G-0 rule', except in case of cycle pair 20,21. That is, this inverse G-0 rule in 
An(T^) and A$(Tm) is violated by about 22 years (a Hale cycle) earlier than the 
corresponding violation of the G-0 rule in Rm- Hence, the studies of variations 
in ^N^m) an d As(T^) may also help for predicting a future even- and odd- 
numbered cycle pair which will violate the G-0 rule, by about 22 years (a Hale 
cycle) advance. Since cycle pair 22,23 satisfies the aforementioned inverse G-0 
rule, hence, probably cycle pair 24,25 will satisfy the G-0 rule, i.e. probably 
cycle 25 will be stronger than cycle 24. 

In Figures 1(a), 1(b) and 2(a) it can be seen that a «44-year periodicity 
exists in Rm also. However, it is not as strong as the corresponding periodicity in 
A-^(T^) and As(T M ) and in their corresponding north-south differences. Rozelot 
(1994) found that in sunspot data only 11-year periodicity is statistically signif- 
icant and all the remaining periodicities are only minor fluctuations. As can be 
seen in Figure 2(a) the epochs of the cycles 15, 19, and 23 are at the minimum 
epochs of the «44-year cycle in <L4(T^), whereas the corresponding epochs of T M 
are at the maximum epochs of the 44-year cycles in SA(T M ). In Figure 6 it can be 
seen that the odd numbered cycles 11, 15 and 19 are followed by the weak even 
numbered cycles 12, 16 and 20. Cycle 7 is also weak and followed by the strong 
cycle 8. However, the cycle pair 8,9 violated the G-0 rule. As pointed out above 
probably the cycle pair 24,25 will not violate the G-0 rule. All these indicate 
that cycle 24 will be weaker than cycle 23. Thus, from the ss44-year periodic 
variations of SA(T^) and SA(T M ) one can infer that the upcoming cycle 24 will 
be weaker than cycle 23. 

From the patterns of SA(T^) and 5A(T M ) we can also infer that the upcoming 
cycles 24 and 28 may be at the beginning (i.e., the end of the current Gleissberg 
cycle) and the ending minimum epochs of the next 44-year cycle, respectively 
(the cycles 32, 36 and so on may be also weak cycles). In both <5A(T^) and 
8A(T M ), and also in Rm, the patterns of the 44-year cycles are considerably 
different from one 44- year cycle to another. In case of 6A(T M ), the pattern of 
its variation seems to be changing in the alternate 80-90-year cycles, indicating 
that the pattern of the next 44-year cycle may be the same as that of the 
current 44-year cycle. Therefore, since the cycle 21 is stronger than cycle 20, we 
can expect that the cycle 25 will be stronger than cycle 24. This is consistent 
with the speculation above, viz. probably the cycle pair 24,25 will satisfy the 
G-0 rule. 
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Figure 6. Plots of annual sunspot number against time. Near the peaks of the cycles 
corresponding Waldmeier solar cycle numbers are shown. 

3. Discussion 

3.1. A Physical Interpretation of the Relationships Found Above 

According to Equations (l)-(4), a large/less value of any of An(T^), As(T M ), 
<5A(T^) and A^g, i.e. a large/less emergence of the magnetic flux (unsigned) 
in the corresponding 0° — 10° latitude interval of the northern or the southern 
hemisphere during Tjjj and T M , of a cycle n implies a large/less emergence of the 
magnetic flux at Tvi of the next cycle n + 1. It should be noted here that the 
epochs T M of cycle n and Tm of cycle n + 1 are around the begin and the end of 
a same polarity half solar magnetic cycle which comprises the declining phase of 
the cycle n and the rising phase of the cycle n+1, whereas the epoch T^ is around 
the middle of the opposite polarity half solar magnetic cycle which comprises 
the declining phase of the cycle n — 1 and the rising phase of the cycle n. There 
exist strong correlations between A^^T^), 8A n (T^) and Ag in (Tj^), whereas 
there is no significant correlation either between j4N.n+i(T,*) and As in (T M ) or 
between 5A n+1 (T£) and 4 s , ra (T M ). This confirms that A N (T£) and SA(T^) lead 
Ag(T^) by about 5 years. It also suggests that the emergence/cancellation of 
the magnetic flux at T^ of cycle n may influence on the emergence/cancellations 
of the magnetic flux at both the T^ of cycle n and Tm of cycle n + 1. The 
emergence / cancellation of the flux at Tm of cycle n may not have an influence 
on the emergence / cancellation of the flux at T^ of cycle n+1. Even it may not 
have an influence on the emergence of the flux at Tm of cycle n+1. That is, 
the high correlation between As,n(Tj^) and Rm,u+i may be due to the influence 
of the emergence/cancellation of the magnetic flux at T^ of cycle n on the 
emergence/cancellation of the flux at both epochs T M of cycle n and Tvi of cycle 
n + 1. The existence of an 'inverse G-0 rule' in A^{T^) and As(T M ) was pointed 
out in Section 2.5. It may be interesting to note here that a similar behavior was 
also found in the latitudinal gradient of the solar rotation determined from the 
sunspot group data (Javaraiah, Bcrtello, and Ulrich, 2005a). Hence, the emer- 
gence/cancellation of the magnetic flux at the epochs T*^ and T-^ of cycle n and 
Tm of cycle n+1, i.e. all the relationships found above, may be related to the long- 
term variations, ss 22 years and longer, in the differential rotation (Javaraiah, 
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2003; Javaraiah, Bcrtcllo, and Ulrich, 2005a; Javaraiah, Bertello, and Ulrich, 
2005b). 

3.2. The Role of the Solar Meridional Flows for the Relationships Found Above 

It is well accepted that the solar dynamo, which seems to be located near the base 
of the Sun's convection zone, generates the solar magnetic field for solar activity 
and the solar cycle (Rosner and Weiss, 1992; Ossendrijver, 2003). The Sun's polar 
fields, solar meridional flows and differential rotation arc important ingredients 
in the dynamo models (Babcock, 1961, Ulrich and Boydcn (2005) and references 
therein). The solar meridional flows can transport magnetic flux and play a 
major role for the magnetic flux cancellation and the polar field reversals (Wang 
(2004) and references therein). As already suggested in Paper I, the relationships 
between A^ tn (T^) and As iTI (7m) an< ^ -^M,n+i may have a physical relationship 
with the solar magnetic cycle and the temporal variations of the solar equatorial 
rotation rate and the meridional flows. Figure 7 shows the variations in the mean 
meridional velocity, v' , of sunspot groups during the odd and the even numbered 
solar cycles determined by superposing the spot group data during cycles 11-23 
according to the years relative to the nearest sunspot minimum epochs 1867, 
1879, 1890, 1902, 1913, 1923, 1934, 1944, 1954,1965, 1976, 1987, and 1997. The 
data reduction and analysis are the same as in Javaraiah and Ulrich (2006). Now 
we have also used the spot group data during cycles 22-23. This increased the 
size of the data considerably. Hence, it enabled us to determine somewhat reliable 
variations in the mean meridional motions of spot groups in the northern and 
the southern hemispheres during the odd cycles (11, 13, 15, 17, 19, 21, 23) and 
the even cycles (12, 14, 16, 18, 20, 22). Tuominen (1952) noticed the existence 
of a difference in the meridional motions of spot groups during the odd and 
the even cycles and suggested the existence of a 22 year periodicity in the solar 
meridional flow. Piddington (1976) assumed that there exist 22-year meridional 
oscillations which change the sign of the angular velocity gradient with respect to 
the relic magnetic field lines causing 22-year solar magnetic oscillations. Hence, 
in order to look for a longer than 11-year periodicity («22 year periodicity) in 
the meridional flow, in Figure 7 the values of the mean meridional motion of 
spot groups during the odd cycles are plotted against the years 1 to 12 and the 
corresponding values of the even cycles are plotted against the years 13-23. A 
vertical dashed line is drawn at the year 12.5 to identify the ends of the odd 
numbered cycles and the beginnings of the even numbered cycles. This figure 
reveals more information than the Figure 2 of Javaraiah and Ulrich (2006), where 
the data from both the odd- and even-numbered cycles are combined. 

As can be seen in Figure 7 there are considerable differences in the mean 
meridional motions of the spot groups in a given hemisphere during the odd and 
the even cycles. For example, at the beginning of an odd cycle the motions are 
equator- ward directions in the 0° — 10° latitude intervals of both the northern 
and the southern hemispheres (the motions in the northern and the southern 
hemisphere are southbound and northbound, respectively). But this is opposite 
at the beginning of an even cycle, i.e. the motions are in pole-ward directions in 
go _ -^go latitude intervals of both the hemispheres. There is an indication that 
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the motion in the 0° — 10° latitude interval of the northern hemisphere rapidly 
changed (in about two years time) from the pole-ward (northbound) direction at 
the end of an even cycle, say n, to the equator-ward (southbound) at the begin 
of the following odd cycle, say n+1. In 0° — 10° latitude interval of the southern 
hemisphere the motion is not significantly different from zero at the end of an 
even cycle n, but it became significantly different from zero and equator-ward 
(northbound) direction at the begin of the following odd cycle, n + 1. These 
opposite directions of the mean meridional motions of the spot groups suggest 
that around the begin of an odd numbered cycle the solar meridional flows 
may transport magnetic flux across the equator causing cancellation of a large 
amount of the magnetic flux around the equator. The equator-ward meridional 
flows may include large contributions of the down flows at the active regions 
in the sunspot latitude belt. The cancellation of the flux is relatively weak or 
even absent around the begin of an even numbered cycle, because the motions 
are pole-ward in 0° — 10° latitude intervals of both the hemispheres. During the 
declining phase of any cycle, until approximately one year before the end of the 
cycle the motions are equator-ward directions in 0° — 10° latitude intervals of 
both the northern and the southern hemispheres. This may cause cancellation of 
some amount of magnetic flux around the equator during the declining phases of 
the cycles. As we have already found in Paper I, the time-intervals Tjjj and Tjjj 
included the epochs when the motion is changed to pole-ward to equator-ward. 

All the aforementioned suggestions of the variations in the mean meridional 
motions of the spot groups are very well consistent with the differences between 
the odd and the even numbered cycles in the values of An,™ K), 5A n (T£J 
and Ag. n (T^j) and the physical interpretations drawn in Section 3.1, above. In 
addition, the overall pattern of the mean meridional motion of the spot groups 
in 0° — 10° latitude interval of the southern hemisphere suggests that about 
9-10 year periodicity is relatively dominant in the solar meridional flow in this 
latitude interval. The corresponding pattern in the 0° — 10° latitude interval of 
the northern hemisphere suggests that a relatively longer periodicity of about 
15-20 years (probably even longer than this, i.e. «30 years), is dominant in the 
solar meridional flow in this latitude interval. 

Here we concentrated mainly on the mean solar cycle variation of the merid- 
ional motions of the spot groups in 0° — 10° latitude intervals of the northern and 
the southern hemispheres. However, as can be seen in Figure 7 the behavior of the 
mean motions in the high latitudes are substantially different from those of the 
low latitudes. It is consistent with the discussion in Paper I on the implications 
of the north-south asymmetries in the latitude distributions of the solar flare 
activity during a solar cycle (e.g., Garcia, 1990), and on the large north-south 
asymmetry in the sunspot activity during the later Maunder minimum (see 
Sokoloff and Nesme-Ribes, 1994). 

It is interesting to note here that during cycle 23 the mean meridional motion 
of the sunspot groups is stronger than that during about last 9-10 cycles. That 
is, during cycle 23 (around maximum epoch) the motion is strongly pole- ward in 
the northern hemisphere and strongly equator-ward in the southern hemisphere 
(see Figures 8 and 9 in Javaraiah and Ulrich, 2006). That is, the overall mean 
motion is northbound and may be responsible for the cancellation of a large 
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Figure 7. Variations in the mean meridional velocity, v' . of the sunspot groups in the northern 
hemisphere (filled circles and solid curves) and the southern hemisphere (open circles and 
dotted curves) during the odd numbered solar cycles (If, f3, f5, f 7, f9, 2f , 23) and the even 
numbered cycles (12, f4, f6, f8, 20, 22). The minima of the odd- and even-numbered cycles 
correspond to year one and thirteen, respectively. Averages are taken over 2-year intervals. The 
vertical dashed line at year 12.5 is drawn to identify the ends of the odd numbered cycles and 
the beginnings of the even numbered cycles. In both the northern and southern hemispheres, 
the positive values of the velocity indicate pole-ward motions and the negative values indicate 
equator- ward motions. The horizontal solid lines represent the zero values of v' . 



amount of magnetic flux during this cycle (mainly during the declining phase) 
and the flux cancellation might be relatively more in the northern hemisphere 
than in southern hemisphere. Hence, the cycle 23 is weak and the activity 
is slightly more in the southern hemisphere than in the northern hemisphere, 
mainly during the declining phase of this cycle. 
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3.3. Comparison of our Prediction with the predictions Based on the Dynamo 
Models 

A large number of forecasting methods (precursors, spectral analysis, non-linear 
dynamics, solar system dynamics, etc.) are used to predict the amplitudes of the 
solar cycles. There have been many predictions for the Rm of cycle 24, but their 
range is very wide (for the list of these predictions see Kane (2007a) and Obridko 
and Shelting (2008)). Therefore, the solar cycle 24 prediction panel found difficult 
to came out consensus on any single value and supports to the two possible peak 
amplitudes, 90 ± 10 and 140 ±20, for the smoothed monthly value of the interna- 
tional sunspot number of this cycle (http://www.sec.noaa.gov/Solarcycle/ 
sc24/statement_01 .HTML). Our prediction is on the low value side. 

Among all the methods of forecasting of the maximum amplitudes of the 
solar cycles, precursor method based on the correlation between the strength of 
a geophysical phenomenon in the declining phase of a cycle and the strength 
of its immediate following cycle, initiated by Ohl (1966), seems to be most 
successful one (Kane, 2007a). Recently, Hathaway and Wilson (2006) and Kane 
(2007a) by using the current strength of the geomagnetic aa index as a precursor 
predicted 160 ± 25 and 142 ± 24, respectively, for Rm of cycle 24. Schatten et 
al. (1978) used for the first time the strength of the polar magnetic field at the 
preceding minimum epoch of a cycle as the precursive indicator for the amplitude 
of the same cycle. This method has a physical background in the sense that 
the strength of the polar field during the minimum of a cycle is an important 
ingredient for dynamo models of the solar cycle. Using this method and the 
current weak strength of the polar field near the end of the ongoing cycle 23, 
measured in Wilcox Observatory, Schatten (2005) predicted 80 ± 30 and Sval- 
gaard, Oliver, and Kamidc (2005) predicted 75 ± 8 for Rm of cycle 24. Recently, 
on the basis of the current strength of the polar field measured in the Mount 
Wilson Observatory, which is weaker than the strength of the polar filed at the 
preceding minimum of cycle 22, Ulrich, R.K. and Boydcn, J.E. also predicted 
that cycle 24 will be weaker than cycle 23 (private communication from Ulrich, 
R.K. and Boyden, J.E., http: //astro. ucla.edu/~ ulrich/). Hathaway and 
Wilson (2004) suggested a large value 145 ± 30 for Rm of cycle 24, using the fast 
equatorial drift rate of sunspot activity (i.e. fast equatorial meridional motion) 
during the declining phase of cycle 22 found by Hathaway et al. (2003). 

Dikpati, de Toma, and Gilman (2006) and Dikpati and Oilman (2006) by 
inputting the data on sunspot group area during the cycles 12-23 to a flux- 
transport dynamo model, predicted a very large amplitude, 150- 180, for Rm of 
cycle 24. One of the main assumptions of their model is a long, about 17-23 
year, magnetic memory (a slow diffusion of magnetic field). On the other hand 
Choudhuri, Chatterjee, and Jiang (2007) and Jiang, Chatterjee, and Choudhuri 
(2007) by inputting the strength of the observed polar field, i.e., by using the 
same data that used by Svalgaard, Oliver, and Kamide (2005), into another 
type of flux-transport model, theoretically confirmed the low values predicted 
by Schatten (2005) and Svalgaard, Oliver, and Kamidc (2005). This model as- 
sumes a short magnetic memory of about 5-year time (a fast diffusion of magnetic 
fields). Cameron and Schusslcr (2007) studied the origin of the predictive skill 
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of the flux-transport dynamo models, theoretically and observationally by us- 
ing both the sunspot area data during cycles f2-23 and the sunspot numbers 
during cycles 1-23, and found that the assumed values for the parameters of 
the aforementioned flux-transport dynamo models arc ineffective, so that the 
predictive power lies in the temporal behavior of the observational data used in 
these models. They have concluded that since the stronger cycles tend to rise 
faster to their maximum activity (Waldmeier effect), the origin of the predictive 
power in the above mentioned flux-transport dynamo models is the shifts of 
the minimum epochs of the cycles due to temporal overlapping of the cycles, 
and hence the predictive skill does not require a magnetic memory. However, 
Dikpati, de Toma, and Oilman (2006) and Dikpati and Oilman (2006) assumed 
a long, 17-23 year, magnetic memory. Hence, a major contribution to their 
prediction for cycle 24 might have also come before the overlapping period of 
the cycles 23 and 24. In addition, recently Dikpati, Gilman, and dc Toma (2008) 
found that the Waldmeier effect is not present in sunspot area. In our case, first 
of all we have used the data on the spot groups in 0° — 10° latitude intervals 
only. Obviously, the epoch Tjjj of a cycle n is quite far away from the overlapping 
period of the cycles n and n+ 1. The T M of a cycle n is just one year away from 
the Tm of cycle n, i.e., it is also considerably far away from the time of the 
overlapping of the cycles n and n + 1. As already found in Paper I, and T M 
of a cycle n are close to the epochs when the polarity reversals of the global 
magnetic fields take place (Makarov et ai, 2001, 2003). Of course, as already 
mentioned above the T M of cycle n and the Tm of cycle n + 1 are in the same 
polarity half magnetic cycle. 

The values we have predicted for Rm of cycle 24, by using Equations (2) and 
(4), are close to the low values predicted by Schatten (2005) and Svalgaard, 
Oliver, and Kamide (2005). However, by using our method the prediction can be 
made by at least four years earlier than their method. Both Dikpati, de Toma, 
and Gilman (2006), Dikpati and Gilman (2006) and ourselves used the same 
spot group data which are taken from David Hathaway's website (we have used 
the raw data files which contain the daily measurements). They have obtained 
a high correlation (r = 0.98) between the simulated and the observed strengths 
of the eight cycles, 16-23. Our prediction is much more statistically significant 
than their prediction, because we have r = 0.98 from eleven data points. Since 
their model needs 17-23 year magnetic memory, hence, in their simulation for a 
cycle n+la major contribution might have come also from cycle n — 1 and even 
from before it, besides a contribution from the cycle n. Their predicted value 
for Rm of cycle 24 is large probably because of a large contribution to it came 
from the large cycle 22 and even also from cycle 21 which is also a large cycle, 
besides a contribution by cycle 23. They have used a constant value, 14 m s _1 , 
for the speed of the meridional flow, whereas the speed of the meridional flow 
may vary considerably (c/., Section 3.2 above and Javaraiah and Ulrich, 2006). 
In addition, they have assumed a profile for cycle 24 whose size is almost equal 
to that of cycle 23. However, they found a large value even without assuming 
any profile for cycle 24 and were able to correctly reproduce the small amplitude 
of cycle 20. Our case is also based on the long-term variations in the sunspot 
activity. But for the prediction of a cycle n + 1 the main contribution comes from 
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around the preceding minimum and near the maximum of cycle n (for predicted 
Rm of cycle 24 the main contribution comes from cycle 23) and that too only 
from the 0° — 10° latitude interval. However, our case does not rule out the 
possibility of some minor contributions of cycle n — 1 and around the minimum 
of cycle n + 1 itself, particularly from the high latitudes. In fact, it may be 
necessary to take into account of such contributions, because the magnetic fields 
at different latitudes during different time-intervals of a previous cycle might 
contribute to the activity at the same or different latitudes during the next cycle. 
If we can include such contributions, then we may get an improvement in the 
corresponding correlation of Equation (4), from 98% to 100%, and a value for Rm 
of cycle 24 which may be considerably different than the values predicted here. 
Recently, Kane (2007b) found a reasonably good correlation (r = 0.89) between 
the sum of the sunspot group numbers in the latitude intervals 10° — 20° and 
20° — 30° of the northern hemisphere during a preceding minimum of cycle n 
and Rm of cycle n + 1, and predicted 129.7±16.3 for Rm of cycle 24. As we 
have already mentioned in Section 2.4, all our predictions are consistent with 
that now the level of activity is at the declining phase of the current Gleissbcrg 
cycle (Javaraiah, Bertello, and Ulrich, 2005b). However, the earlier result itself 
is yet to be confirmed. 

3.4. Comparison of our Prediction with the Predictions Based on the Spectral 
Analyses and Magnetic Oscillations Models 

The properties of the solar cycles can also be explained on the basis of superpo- 
sitions of the dominant modes of the Sun's global magnetic oscillations whose 
frequencies equal to that of the solar magnetic cycle and a few harmonics of 
it (e.g., Bracewell, 1988; Stcnflo, 1988; Gokhale and Javaraiah, 1990; Gokhalc 
and Javaraiah, 1995; Juckctt, 2003). Such a model has a strong predictive power 
(e.g., Gokhale and Javaraiah, 1995). A number of authors used various spec- 
troscopic methods and predicted around 100 for Rm of cycle 24 (e.g., Kane, 
1999; Echer et ai, 2004). Recently, Hiremath (2006) modeled the solar cycle 
as a forced and damped harmonic oscillator. Hiremath (2008) by inputting the 
amplitudes, frequencies, and phases of 22 cycles (1755-1996) derived from the 
aforementioned model into an autoregressive model, predicted the periods and 
the amplitudes of the next fifteen solar cycles. His prediction for Rm of cycle 24 is 
110 ± 11. It may be interesting to note here that some cycles, including cycle 23, 
contain double peaks. Gnevyshev (1967) interpreted this property of the solar 
cycles as an 11 year cycle consisting of two processes (waves ?) with different 
physical properties and the shape of the cycle depends on the way these processes 
overlap. That is, the occurrence of a double peak cycle may be mainly related to 
the difference between the phases of the processes (waves) involved. That is, a 
large difference in the initial phases of the waves may be responsible for a small 
and broad or double peak, and also long period, cycle. A number of authors also 
explored the possibility that the solar cycle is a sum of two periodic functions 
(e.g., Merzlyakov, 1977). The variations in SA(T^) and A$(T M ) may represent 
the two dominant magnetic waves (quadrupole and dipole components of the 
global magnetic oscillations) in the Sun, whose superposition may be responsible 
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for the long-term variations in the amplitudes of solar cycles. The A^iT^) and 
5A(T^) lead As(T M ) by about 4 years. The difference between Tj£ of cycle n 
and Tm of cycle n + 1 is 14-18 years, and that between T M of cycle n and Tm of 
cycle 7? + l is 9-12 years, which may represent the periods of the aforementioned 
waves. Such periodicities are found to exist also in the solar differential rotation 
determined from sunspot group data (Javaraiah and Gokhale, 1995; Javaraiah 
and Komm, 1999; Javaraiah, 2003; Javaraiah, 2005; Brajsa, Ruzdjak, and Wohl, 

2006) . The existence of ~ 22-year periodicity is detected also in the solar rotation 
determined from the data on large-scale magnetic fields measured using the 
Ha line, magnetographic observations, and spectral-corona observations (Tlatov, 

2007) . We think the combined effect of the Sun's rotation and the inclination of 
the Sun's equator to the ecliptic may have a major role in all the relationships 
found above. 



4. Summary and Conclusions 

Using the values of the sums of the areas, An and Ag, of the spot groups in 
go _ latitude intervals of the northern and the southern hemispheres and 
during the time intervals Tjjj and T M of cycle 12-23 determined in Paper I, we 
have found the following: 

1. A^(T^) varies strongly with a period of about 44 years (double Hale cycle), 
whereas there exits a weak 80-90 year Gleissberg cycle in As(T£). 

2. A$(T M ) is also having a 44-year periodicity, whereas An(T m ) has a somewhat 
different periodicity of about 33-44 years. 

3. There exist statistically significant north-south asymmetries in the aforemen- 
tioned sums of the areas of spot groups, i.e., there exist statistically significant 
differences <L4(T*) = A N (T£) - A S (T*) and SA(T M ) = A N (T M ) - A S (T M ). 

4. 5A(T*A and SA(T M ) vary approximately in opposite phase. Both these vary 
strongly with a periodicity of about 44 years. 

5. The patterns of the 44-year cyclic variations in SA(T^) and 5A(T M ) strongly 
indicate that the upcoming cycle 24 will be weaker than cycle 23 and it will 
be the beginning minimum epoch of the next 44-year cycle (it may be also 
the end of the current Gleissberg cycle). The next such a minimum may be 
represented by cycle 28. 

6. As in the case of A N (T^) and As(T M ), the north-south difference, SA(T^) of 
a cycle, correlates well with Rm of its immediate following cycle. By using 
this relationship it is possible to predict the amplitude of a cycle by about 
13-year advance, with an improved accuracy. Using this relationship we have 
obtained 103 ± 10 for the amplitude of the upcoming solar cycle 24. 

7. The correlation between the sum SA(T^) + As(T M ) of a cycle and Rm of its 
immediate following cycle is very high. By using this relationship it is possible 
to predict the amplitude of a cycle by about 9 years in advance with a high 
accuracy. Using this relationship we have obtained 87 ± 7 for Rm of cycle 24, 
which is about 28% less than the Rm of cycle 23. 
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8. In the cycle-to-cycle variations of A-^(T^) and As(T^) there is an indication 
of the upcoming cycle pair 24,25 will not violate the G-0 rule, i.e. it seems 
cycle 25 will be stronger than cycle 24. 

9. There exists a considerable difference between the variations in the mean 
meridional motions of the spot groups during the odd numbered and the 
even numbered cycle. These variations suggest that there exist about 9-10 
year and 15-20 year periodicities in the mean meridional motions of the 
spot groups in 0° — 10° latitude intervals of the southern and the northern 
hemispheres, respectively. 

10. The spacial (latitudinal) patterns of the aforementioned variations in the 
mean meridional motions of the sunspot groups suggest that the solar merid- 
ional flow may transport magnetic flux across the solar equator and po- 
tentially responsible for all the relationships found here. Consequently, our 
results are not only highly statistically significant and useful for an accurate 
prediction of the amplitude of the upcoming cycle 24, but also have many 
implications for understanding the underlying physical process responsible 
for the solar cycle. 

Acknowledgements I thank the anonymous referee for useful comments and 
suggestions. I also thank Professor Roger K. Ulrich for fruitful discussion, and 
Dr. R. W. Komm, Dr. D. H. Hathaway and Dr. K. M. Hiremath for useful 
comments. I acknowledge the funding by NSF grant AM-0236682. 



References 

Antonucci, E., Hoeksema, J.T., Scherrcr, P.H.: 1990, Astrophys. J. 360, 296. 
Atac, T., Ozgiic, A.: 1996, Solar Phys. 166, 201. 
Babcock, H.W.: 1961, Astrophys. J. 133, 572. 

Ballester, J.L., Oliver, R., Carbonell, M.: 2005, Astron. Astrophys. 431, L5. 

Bracewell, R.N.: 1988, Mon. Not. Roy. Astron. Soc. 230, 535. 

Brajsa, R., Ruzdjak, D., Wohl, H.: 2006, Solar Phys. 237, 365. 

Cameron, R., Schiissler, M.: 2007, Astrophys. J. 659, 801. 

Carbonell, M., Oliver, R., Ballester, J. L.: 1993, Astron. Astrophys. 274, 497. 

Chang, H.Y.: 2008, New Astron. 13, 195. 

Choudhuri, A.R., Chatterjee, P., Jiang, J.: 2007, Phys. Rev. Lett. 98, 131103. 
Dikpati, M., Gilman, P. A.: 2006, Astrophys. J. 649, 498. 

Dikpati, M., de Toma, C, Gilman, P.A.: 2006, Geophys. Res. Lett. 33, L05102. 
Dikpati, M., Gilman, P.A., and de Toma, G.: 2008, Astrophys. J. 673, L99. 
Duchlev, P.I., Dermendjiev, V.N.: 1996, Solar Phys. 168, 205. 

Echer, E., Rigozo, N.R., Nordemann, D.J.R., Vieira, L.E.A.:2004, Ann. Geophys. 22, 2239. 
Fairbridge, R.W., Hillaire-Marcel, C: 1977, Nature 268, 413. 
Garcia, H.A.: 1990, Solar Phys. 127, 185. 
Georgieva, K. : 2002, Phys. Chem. Earth 27, 433. 

Georgieva, k., Kirov, B. : 2003, In: Wilson, A. (ed.), SOHO 12/GONG+ 2002, Local and 

Global Helioseismology: The Present and Future, ESA SP-517, 275. 
Georgieva, K., Kirov, B., Javaraiah, J., Krasteva, R.: 2005, Planet. Space Set. 53, 197. 
Gigolashvili, M.Sh., Japaridze, D.R., Kukhianidze, V.J.: 2005, Solar Phys. 231, 23. 
Gnevyshev, M. N.: 1967, Solar Phys. 1, 107. 
Gnevyshev, M.N., Ohl, A.I.: 1948, Astron. Zh. 25, 18. 

Gokhale, M.H., Javaraiah, J. : 1990, Mon. Not. Roy. Astron. Soc. 243, 241. 



javaraiah_msrev.tex; 8/03/2009; 7:44; p. 22 



PREDICTION OF THE AMPLITUDE OF A SOLAR CYCLE 



Gokhale, M.H., Javaraiah, J. : 1995, Solar Phys. 156, 157. 

Hathaway, D.H., Wilson, R.M.: 1990, Astrophys. J. 357, 271. 

Hathaway, D.H., Wilson, R.M.: 2004, Solar Phys. 224, 5. 

Hathaway, D.H., Wilson, R.M.: 2006, Geophys. Res. Lett. 33, L18101. 

Hathaway, D.H., Wilson, R.M., Reichmann, E.J.: 1999, J. Geophys. Res. 104, 22375. 

Hathaway, D.H., Nandy. D., Wilson, R.M., Reichmann, E.J.: 2003, Astrophys. J. 589, 665 

(Erratum: 602, 343 (2004)). 
Hiremath, K.M. : 2006, Astron. Astrophys. 452, 591. 
Hiremath, K.M. : 2008, Astrophys. Space Sci. 314, 45. 
Javaraiah, J.: 2003, Solar Phys. 212, 23. 
Javaraiah, J.: 2005, Mon. Not. Roy. Astron. Soc. 362, 579. 
Javaraiah, J.: 2007, Mon. Not. Roy. Astron. Soc., 377, L34 (Paper I). 
Javaraiah, J., Gokhale, M.H.: 1995, Solar Phys. 158, 173. 
Javaraiah, J., Gokhale, M.H.: 1997, Solar Phys. 170, 389. 
Javaraiah, J., Komm, R.W.: 1999, Solar Phys. 184, 41. 
Javaraiah, J., Ulrich, R.: 2006, Solar Phys. 237, 225. 
Javaraiah, J., Bertello, L., Ulrich, R.: 2005a, Astrophys. J. 626, 579. 
Javaraiah, J., Bertello, L., Ulrich, R.: 2005b, Solar Phys. 232, 25. 

Jiang, J., Chatterjee, p., Choudhuri, A.R.: 2007, Mon. Not. Roy. Astron. Soc. 381, 1527. 

Joshi, B., Pant, P.: 2005, Astrono. Astrophys. 431,359. 

Juckett, D.: 2003, Astron. Astrophys. 399 , 731. 

Kane, R.P.: 1999, Solar Phys. 189, 217. 

Kane, R.P.: 2007a, Solar Phys. 243, 205. 

Kane, R.P.: 2007b, Solar Phys. 246, 471. 

Knaack, R., Stenflo, J.O., Berdyugina, S.V.: 2005, Astron. Astrophys. 438, 1067. 
Li, K.J., Yun, H.S., Gu, X.M.: 2001, Astron. Astrophys. 368, 285. 

Li, K.J., Wang,. J.X., Xiong, S.Y., Liang, H.F., Yun, H.S., Gu, X.M.: 2002, Astron. Astrophys. 
383, 642. 

Makarov, V.I., Tlatov, A.G., Sivaraman, K.R.: 2001, Solar Phys. 202, 11 

Makarov, V.I., Tlatov, A.G., Sivaraman, K.R.: 2003, Solar Phys. 214, 41 

Merzlyakov, V.I.: 1997, Solar Phys. 170, 425. 

Obridko, V.N., Shelting, B.D.: 2008, Solar Phys. 248, 191. 

Ohl, A.I.: 1966, Soln. Dannye 12, 84. 

Osscndrijvcr, M.: 2003, Astron. Astrophys. 11, 287. 

Piddington, J.H.: 1976, In: Bumba, V., Kleczek, J. (eds.), Basic Mechanisms of Solar Activity, 
IAU Syrap. 71, 389. 

Rosner, R., Weiss, N.O.: 1992, In: Harvey, K. L. (ed.), The Solar Cycle, ASP Con}. Ser. 27, 
511. 

Roy, J.R.: 1977, Solar Phys. 52, 53. 
Rozelot, J. P.: Solar Phys. 149, 149. 
Schatten, K.: 2005, Geophys. Res. Lett. 32, L021106. 

Schatten, K.H., Sherrcr, P.H., Svalgaard, L., Wilcox, J.M. : 1978, Geophys. Res. Lett. 5, 411. 
SokolofT, D., Nesme-Ribes, E.: 1994, Astron. Astrophys. 288 293. 
Stenflo, J.O.: 1988, Astrophys. Space Sci. 144, 321. 

Svalgaard, L., Oliver, E.W., Kamide, Y.: 2005, Geophys. Res. Lett. 32, 021664. 
Swinson, D.B., Koyama, H., Saito, T.: 1986, Solar Phys. 106, 35. 

Temmer, M., Rybak, J., Bendik, P., Veronig, A., Vogler, F., Otruba, W., Potzi, W., Hanslmeier, 

A.: 2006, Astron. Astrophys, 375, 1049. 
Tlatov, A. G.: 2007, Astron. Lett. 33, 771. 
Tuominen, J.: 1952, Z. Astrophys. 30, 261. 
Ulrich, R.K., Boyden, J.E.: 2005, Astrophys. J. 620, L123. 
Verma, V.K.:1993, Astrophys. J. 403 , 797. 
Wang, Y.-M. : 2004, Solar Phys. 224, 21. 
Yi, W.: 1992, J. Roy. Astron. Soc. Canada 86, 89. 

Zaatri, A., Komm, R., Gonzalez, H.I., Howe, R., Corbard, T.: 2006, Solar Phys., 236. 227. 



javaraiah_msrev.tex; 8/03/2009; 7:44; p. 23 



javaraiah_msrev.tex; 8/03/2009; 7:44; p. 24 



